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carborundum when used with a “whisker” may either 
decrease or increase with the increase of temperature of 
the crystal. We have observed these features, and it 
seems likely that Sen (Ind. Jour. Physics, vol. x. 1936) 
has also noticed them. 
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XIII. The Diffraction of Electromagnetic Waves from an 
Electrical Point Source round a Finitely Conducting 
Sphere, with Applications to Radiotelegraphy and the 
Theory of the Rainbow.—Part I. By BaLTH. VAN DER 
Por and H. BREMMER, Natuurkundig Laboratorium der 
N.V. Philips’ Gloeilampenfabrieken, Eindhoven, Holland*. 


1. Introduction. 


T paper concerns the calculation of the electro- 

magnetic field due to an oscillating electrical dipole 
placed outside a finitely conducting sphere. This is 
one of the fundamental problems of radiotelegraphy 
and has aroused new interest by the coming of television. 
Here the radiating dipole represents the transmitter 
antenna, and the absorbing sphere the earth. In our 
calculations the atmosphere is taken homogeneous and 
any effects due to the ionosphere are therefore neglected. 
(The chief contents of Part I. of this investigation have 
been communicated by Balth. van der Pol in a lecture 
before the Institute of Radio Engineers in New York in 
November 1936.) 

The rigid solution for a given wave-length, based on the 
Maxwellian equations with proper boundary conditions, 
leads to an infinite series of spherical harmonics with 
coefficients containing twelve Bessel functions, which 
has been known already for a long time and which is given 
below as formula (1). 

As in all diffraction problems the character of the solu- 
tion is determined by the ratio of the dimensions of the 
obstacle to the wave-length A, in case 27a/A, where a is the 
radius of the sphere. When this ratio is small, as, for 
instance, in the problem of the scattering of light by sub- 
microscopic particles, the series (1) converges so rapidly 


* Communicated by Dr. Balth. van der Pol. 
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that a few terms only give a sufficient numerical approxi- 
mation. For instance, this approximation led Lord 
Rayleigh ® to discover his A~* law for scattered light, and 
it also is at the base of Debye’s investigation of the light 
pressure ), 

In the radio application of our problem this ratio is 
of the order 10° to 108, while the corresponding ratio 
in the theory of the rainbow is determined by the dimen- 
sions of the raindrops and the wave-length of light, leading 
to a value of 27a/A of the order of 104. Under these 
circumstances the series (1) converges in general extremely 
slowly, the main contribution being given by the terms 
where n is of the order 2ia/A. Lord Rayleigh @ evaluated 
for an analogous acoustical problem with a similar, but 
slightly different boundary condition, the field at the 
surface of the sphere, the point source being also on the 
surface, by simply adding numerically the important 
terms for the cases 27a/A=2 and 10. Later on, A. E. H. 
Love @) even undertook the laborious task of evaluating 
the series (1) for the same case of radiotelegraphy with 
2ra/X=8000, by numerical addition of the terms of the 
series taken in groups. In his paper an historical summary 
of the problem up to 1915 is also given. 

Considerable progress in the present problem was 
subsequently obtamed by G. N. Watson ©), who trans- 
formed (1), with the aid of a contour integration, into 
another much more rapidly converging series (15), which 
was further considered in more detail with a view to the 
application to radio by one of us‘®, Also Laporte’s 
investigation uses this series, and a clear exposition by 
A. Sommerfeld of Watson’s method is to be found in 
Frank-Von Mises, Differential Gleichungen der Physik, 
Ii, 2. vermehrte Auflage, pp. 964-977. Finally, T. L. 
Eckersley “) treated the problem with his method of 
phase integrals, while P. S. Epstein ® simplified the 
treatment by an application of Huygens’s principle. 

The present investigation attacks the problem in two 
quite different ways :— 


I. The first one, extending Watson’s treatment, is 
more appropriate for investigating the disturbance 
beyond the optical horizon of the point source and is 
pushed to considerable numerical details while considering 
.a.0. in the radio case the effect of the height of the trans- 
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mitter above the earth. E.g., we find that for practical 
cases of modern television transmission (A=7 meter) the 
attenuation of the radiowaves due to the finite conduc- 
tivity of the earth dominates by far the attenuation due 
to diffraction. Further, when moving along the surface 
of the earth, the optical horizon is found not to be a very 
critical point as regards the amplitude of waves of the 
order of A=7 meter, more or less sharp shadows only 
occurring for much shorter wave-lengths. Also it is 
found that for this wave-length the factor by which the 
field on the surface of the earth and beyond the horizon 
is increased due to raising the transmitter to a certain 
height is, as a first approximation, independent of the 
distance. Moreover, the numerical results show that, 
as in the analogous problem for a plane boundary, the 
conductivity of the earth has a very pronounced influence 
on the attenuation of the waves. 


II. The second method also starts with the rigorous 
harmonic series (1), being a solution of the wave equation. 
It appeared possible to split up this series uniquely in 
a series of an infinite number of other series, each of which 
may be considered to represent the rigorous expression 
of what can be visualized as rays undergoing successive 
reflexions at the inner surface of the sphere. An approxi- 
mation of each of these series clearly shows their ray- 
character in complete accordance with the classical 
Airy ray-treatment of the theory of the rainbow. The 
present investigation therefore provides the link between 
the exact wave solution and the approximate ray solution 
of the rainbow extending the latter to the cases of a 
light source and observer at finite distances. 


2. The Development in Spherical Harmonics. 


Using spherical coordinates, let (fig. 1) the electrical 
dipole be situated at Q(b, o), the observer at P(r, 3), the 
radius of the sphere being a. The analytical problem 
consists of finding the function JT, representing the ampli- 
tude of the Hertzian vector divided by r, having a radial 
component only, and satisfying the following differential 
equations and boundary conditions (where k? , has the 
usual meaning) : 

(1) (4+47=0 (r>a), 

(4--k?)7=0 (r<a), 
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(2) 2 (rf) and kIT continuous at r=a, 
r 


ik1R 
(3) a singularity at Q(b, o) of the form TR which is 

t 
further called. the “ primary field.” This means that in 
the outer space H =H, preo where the “secondary 


Fig. 1. 


Coordinates used. 


field” hee has no singularities. The electric and 
magnetic fields are given by 
a 1 @ 
BE E 28 - : ZE : 
E= (K+ ga) 0D; Ey= = ars, (rll): 
E,=0, H,=Hy=0; Hy=— Sa 2; 


w 


g €i, 20*-+ 10}, 9 
13 c2 


The known solution in the form of harmonic series is : 


a g% 
Tye =I y+ see — ik R 
E 2n-yR, AeA) cor(e.b)e(kyr)P,(0089) (ra), 
a=0 ‘Se (ka) 
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Tor k E (20+ 1)(1+-By) 22D tk Bin egr)Pn(co8 9) 
n=0 


key palka) 
(r<a), . . (1b) 
where the symbols have the following meaning : 
ps ld y ee 
Dy) — T HO =| — Z ae 
on (x)= Jz. HY ,@) x x dz ir > 


Pasa / THe, (wma (—) ZY (HE) 
Pn(%)=$ {CO e) +62 (x)} 
=V gme (a (2), 
(H.=Hankel function), and where 
gles] + at CZOIR 


E pese) [zz later] 


& 

- (2) 
It may be remarked that (1a) is symmetrical inr and b, 

expressing the well-known property of reciprocity, 

stating that transmitter and receiver may be interchanged. 


x 


3. Transformation of (1) to a Sum of a Series of 
Residues and an Integral. 


For the derivation of (1) it is necessary to develop 
first the primary field in a harmonic series. It is well 
known that two different series are obtained according 
to whether r>6 or r<b, viz. : 


elak 


CT = E (MLD (kyryq(k,)P,(c08 9) (rb), (3a) 
ik,R n= 


tkiR =o 

one X (Qn +1) (k,b Wn (Kyr)P,(cos 8) (r<b). (3b) 

tk,R n=0 

If we express (la) completely as a harmonic series 
by substitution of (3b) into (1@) and considering the 
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total field at the surface of the sphere only {r =a), while 
making use of the Wronksian : 


PAPAL 0. 


we obtain 
((k,b) 1 


T= yey 2 (nD Fae a) N, 


P, (cos 3), . (5) 
where 


Nate, y= [1H log {atoa} ] 


r= kya 


ld 1 
—j|-—lo a . (6 
[pay ee ant |) 
Following Watson, we transform (5) into a continuous 
integral over n leading to 


1 ndn 4 0(kyb) 1 
Tia= wath cos (n7) S 1 AEA Na- P,,—10{c08 (a=) 


yal i> ndn Kaetah ( S de ) 
(ka)? o Cos (n7) gu. 12 (kya) Naiz Naga 12 
x P,-12{e08 (7—S)}.. (7) 


Here the contour of integration L in the complex 
n-plane lies wholly in the first quadrant and encloses all 
poles lying in that quadrant {see fig. 2c), This path of 
integration actually makes the sum of the two integrals 
(7) identical with the series (5), which may be clear by 
the following transformation:—First the series (5) is 
transformed into an integral with the path of integration 
formed by small circles round the poles n =1/2, 3/2, 5/2 . . 
of the function 1/cos (n 7) (see fig. 2 a, where ng, ny, na... 
are the poles of the integrand apart from the factor 
1/cos (nz)). Further, this path can be transformed to the 
loop A of fig. 2 b, extending from œ —ietoco-+ie. As the 
integral becomes exponentially zero when || tends 
to oo, the integral with integration path A is the differ- 
ence of the contour integral B of fig. 26 and an integral 
along the imaginary axis. Further, it appears that the 
integrand has no poles in the minus. first quadrant, so 
that B may be contracted to the contour L of fig. 2c, 
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whereas the integral along the imaginary axis may be 
transformed to the second integral of (7) along the positive 
imaginary axis. 
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T k -b 


Fig 2a 4 
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Fig 2b 
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Fig 2c 
Transformation of the path of integration in the n-plane, 


Apparently the second term of (7) vanishes when 
Ny—12 is an even function of n. Now the first term of 
(6) is even in n, whereas the second term, for the case of 
radiotelegraphy, is to a high approximation equal to 


d 
igote] 
L2 
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which is also an even function of n ; this approximation 
is a consequence of the fact that in practice | e2™z f 
(corresponding to the attenuation of plane waves which 
have travelled over a distance equal to the diameter of the 
earth) is very small, due to the great absorption in the 
earth. Therefore, as a first approximation, we may 
consider N,_15 as an even function of n and omit the 
second term of (7). Watson (l.c. p. 99) showed 
also for a numerical case that this term is negligible 
compared with the first term, and further on we shall 
neglect this term. 

The first term of (7) enables us to obtain a very physical 
interpretation of a typical possibility inherent to our 
spherical problem and which has also been noticed by 
Watson. For if we develop 1/cos (n7) as 


-l gem" g ginmensy, (gy 
cos (n7) m=0 


valid over the whole integration path L, we can write 
(T) as 


Iot E Tm ee a eee a (9): 
where 
i 2em n dn etmn(Qm+1) Eag (kb) T z 
(ka)? L gD 2 (ka) N. —12 


x P,,-1p{cos (w—4)} Be ii (10} 


We are now able to represent this integral as the sum 
of the residues of all the enclosed poles being the zeros 
n, Of N,_1e only, and therefore may write it as 


dgicinmes (0. (kb) 1 
f= E nemm) ee -e 
™ (kaf sa LO 5 (kya) (z) 

on n=, 


x Pn,-1{ cos (7—S)}. . (11) 

Further, for all values of 3 not too near 0 or r, and, 

as the modulus of all poles n, is of the order ka (as will 

be shown later on), the well-known asymptotic formula 

for the spherical harmonics of complex order with positive 
imaginary part 


Pa- {cos (r—S)}~ —e-ine—S)+imit, (12), 


1 
vrn sin% 


Electromagnetic Waves from an Electrical Point Source. 149 


may now be introduced as an approximation in the still 
rigorous expression (10), yielding 


2V Fr ertaro 


ge ee ee 
(kav sins 


s=% one —12(e,b) o ng? ein(St2mm), (13) 
s=0 one —12( Key) ( ON, -1 2) 
on R=Ns 


The above mentioned physical interpretation is given 
directly by the last factor of (13), viz. : 


eitg(9+ 20m) — elias— B94 2m) . . . 04 
putting 
Ng =, Hib; (8>0). 


Hence (13) represents a wave travelling round the 
sphere with a phase velocity and attenuation determined 
by æ, and ß, respectively. But the presence in the 
exponent of the combination $+2rm clearly shows 
that the term H, of (9) is representative of a wave 
which has travelled from the source to the receiver only 
after having made m complete revolutions round the 
sphere. This explanation becomes especially vivid for 
the case of radiotelegraphy, where the transmitter is near 
the surface of the sphere. It will be clear that the con- 
tribution to the field at the surface of the earth by a 
wave which has travelled several times round the earth, 
is extremely small compared with the wave coming 
directly from the transmitter to the receiver (m—O). 
In Watson’s analysis the first term IJ, of the series (9) 
is retained only, and we will follow the same procedure in 
Part I. of this paper, so that we will further consider the 
approximation 


2v Pred tin 
Whit] hy 
ee (ka)? sin § 


(1) 
see e cekib) "gins, (15) 
smo EO 12 (kua) (Bs) 
on 


R= hy 
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4. Approximations for the Bessel Functions. 


It will be clear that for numerical interpretation of (15).,. 
a knowledge of the poles n, and of the residues is necessary. 
We will therefore make use of two approximations : 


(a) The Debye approximations, as given in full detail 
in the second edition of Jahnke-Emde, Funktionentafeln, 
pp. 204-207. These become in our case, where n/z and z 
lie in the first quadrant : 


gy (z)= are ie Oy HE o Lye (2) 
et 2 — 
: == 
inl (n f ł)z V, are ae 
x cosd 7 +(n+4 ye vari X ». (16a) 
ioe) = Jz HE vol) 
1 


OM wren ee 


i {n+ 1/2) af 


inj/d-a +12) ui y ae 
E 


(16b) 

where the root a ai „—1 is defined as lying in 
(n+4) 

the same quadrant as z/(n+4). 


(6) The approximations reducing the Bessel functions 
of any order to these functions of order 1/3. For the 
orders of magnitude obtaining in our problem, these are 


Je 
mt 7 sinjsin |" (ue 13.2 
m 
aBC w 


2 
-r2 < arga / Z —1 <n/2. 


where 
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The formule (16) are suitable everywhere except for 
z near n, while (17), valid only for n>>1, is of special 
use when z is near n. It is further easy to verify that 
the first term of the asymptotic development of Hy); 
occurring in (17) leads to (16). 

Moreover (16a) clearly shows the position of the 
zeros of {® (z) at constant n, these being determined by 
the condition that the cosine becomes zero, similar roots 
not occurring in the quadrant considered in the case of 
[2 (z). Again, these zeros of (@ (z) separate the two 
regions where one or the other of the exponential parts 
of the cosine function dominates, these zeros requiring 
obviously 


p “ziin 1/2) ras: 
Ing ri +(n-+4){ visi b =0, . (18) 


ul 


where I denotes the imaginary part of its argument. 

The close relationship between (16) and the well-known 
W.K.B. approximation is at once apparent, remembering 
that zZ,(z) satisfies the differential equation 


Detd MED y ato. . . (19) 


5. The Limiting Case of a Plane Boundary (a> ). 


In this paragraph we will show how the well-known 
solution of the similar problem, but for a plane boundary, 
can be obtained as a limiting case of the spherical problem. 
We therefore start with (1 a) in the form 


oe, Palkia) y a) 
HMes= 2. (2n+1)R, 7 (ka) CO(k,b)CM(kyr)P, (cos 5) 


(r>a), . (la) 


and have to consider the limit of this expression for 
a> while keeping constant the three quantities 


b—a =z, (height of source above boundary), 
2—a =z, (height of receiver above boundary), 


aS =p (horizontal distance between source and receiver). 
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Using the method of §3 in order to transform (la) into 
an integral, one finds 


F d 
w lho 
B 


cos (nr) 


(k) p  afeos(m—3)}. (20) 


Now putting n =@à, = N the limit for aco , while 
keeping A constant, one finds, for the different factors 
of the integrand of (19), 
(a) 
24/2 E N E 
Lim Ra-1,2= kı A ky a x ae . (21) 
nae, k ev 2 Ttk 24/2 kes 27 

Here we have used the limits of the logarithmic derivates 
of (17) in the form 


d D 1 nm £2 3/2 
= jj ae Off = 
dz log fe? (@)} 2z Fr z (r 1) 


+ gethin Z 


n nè Hed n ey a2)? 


7 2 T 
=s <argy/ 7 ~1<2. 
In our limit z/n has the constant value k/A, but the 


argument of Hı; and Hy, is tending to +æ ,and we 
obtain 


where 


fan = $ tog fet! (2)} = F pV A2, 


th 


where V A?—k? has always a positive real part. Further 
we have 


Lim g,(ka)=} Lim £2(ka), 


because k lies in the first quadrant and therefore the 
part vanishes, so that 


Lim z log{ahl2)}=4Lim # log (202). 
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The limit of R,,_ 1/2. then follows directly from the definition 
(2) of R,. 

(b) 
1 


Lim © (k {ka= — : 
se ra yal 1), 12 18) ikay k 


because, as in (a), we first have 
tam as 12l ka)fn—i 2( kya) =} Lim gr 1 p(ka) rA alkia), 
a> o a> 
and then can apply (16), using in (16 a) only the dominant 
term of the cos-function. 
(c) 
Lim g (kyb) _ 
teasel a) 
This limit follows directly from the dominant term in 
(16 a), because 
E yolkb) Y a vi a— n 
ED elka) Vo 7k 2b? —n2” 
and this expression leads at once for a>% to the above 
limit. 
(d) Similarly we have 


— Vat kg, 


, (AL beg apk A — du 
ia {pon garki 4/ue—1 A 
> 


Lim CO (kar) ae AT Eizy 
eta 
(e) 
Lim P„_12{c08 (r—3)}_ Lim P,_12{cos (m —3)} 
In)>0 cos (nz) lino cos (nz) 


where I(n) means the imaginary part of n. 
We may prove this by using the relation 
P,,_12 {008 (m —5)} 
cos (17) 


=1g(27)Py,1,2(08 8) + = Qax (0084), 


which is valid for unrestricted n but real $. For n=aa 
tending to © the functions P,,_12 (cos $) and Q,_1,2 (cos 5) 
have each the same limit for I(n)>0 and I(n)<0, but 
for tg (nz) we have 

Lim tg(nm)=-iresp. for I(n) Z 0. 


|n|oo 
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Therefore 
Lim P,,_;o{cos(7--S)}_ Lim P,,_, o{cos (w—*)} 
Ie) 0 cos (nz) limo cos (nr) 
—2i Lim P,_19 (cos $) 


how 


=2i L'm Ppap (cos =) ==2tJ (àp). 


n> w 


After substitution of the limits a—d in (20), which are 
valid both for the part of B above and the part below the 
real axis, we can reduce the path of integration to the real 
axis after applying (e), and finally obtain 


1 i MA (keV Eh-he V Xk?) 


T= iki 


xe- Viinit) J (Ap). 


This coincides with the expression given as (1) in a former 
paper °°), where the primary field was defined as e8/R 
instead of e¥/ik,R as here. Hence a complete agree- 
ment is obtained, including a finite height of both point 
source and receiver and a finite conductivity of the two 
media. 


6. Development in Characteristic Functions. 


It may be noted that the solution used so far is quite 
distinct from the well-known method of development 
in Characteristic Functions, as, e. g., is a standard method 
for string problems and in wave mechanics. 

In our problem we have two differential equations for 
r>a and r<a; before applying the theory of Character- 
istic Functions, we shall first consider the more general 
equation of the Sturm-Liouville type 


{A4+g(a, y, 2) fT =f(a, y, 2). 


This equation has Characteristic Functions defined by 
the following properties : 


1. (A+A;*9)$;=0, 
2. ¢; finite for r=0 and r=% . 


The expansion of the solution in Characteristic Functions 
as described for instance in chapter v., §2 of Courant und 
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Hilbert, Methoden der Mathematischen Physik, i., here 
becomes 


mP)= a ietdr apy, ee (23) 


Now let g(x, y,z) degenerate to the discontinuous 
function g=k,? for r>a and g=k for r<a, gI and 
lər (rII) remaining continuous for r=a; then the new 
Characteristic Functions are determined by the conditions 

(1) (4+4,7A7)¢;=0 (r>a) j 
(A-k p= (r <a), 


(2) kd; and 3 „("$;) continuous for r=a,}. (4) 


(3) me] TA pdh A Maas: 


We now proceed by a ieod which has been fully 
expounded by A. Sommerfeld “Y in a classical paper, 
where the primary field is cos (k,R)/(ik,R), belonging to 
a point source of stationary waves with a singularity at Q ; 
this field Jeads to a three-dimensional impulsive function 
4718(R)/k, at Q for f(x, y,z). Hence substitution in (23) 

ields 


yiel 
i ár n PQP). 
== 2; 3 e] 
Mios ik, (1—28) 
As the Characteristic Functions should everywhere be 
finite and single-valued, they are of the form 


$ — Ag, EPA, jkır)Pn (cos $) (r>a), ) 
FU Bp sthalAn, kar) Pa (cos) (r<a), J” 


where » is a positive integer. 
The boundary conditions further require 


kP An, EP On, kaa) =k By, wh, (An, hod), - (27a) 
d 7 d 
Ans (a EPD] pa SE Lae] 


5 tele, (25) 


(26) 


3 
Z=Ay, jka 


- - (275) 
so that 
kParLPy, skye) apal Ay, seat) 
d at 
ere]. [Stone al, 


or N,(An, jka, ht) =0, (28) 
(see (6)). 
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The difference between the method used. so far and the 
present method of Characteristic Functions now becomes 
very clear by the fact that, although both methods consider 
the zeros of N,—which is a function of n and the ka’s— 
the first method considers the zeros as a function of 
n for constant ka’s, whereas the method of Characteristic 
Functions treats the zeros as a function of the ka’s, 
while keeping the integer n constant. 

Introducing the normalization condition (3) of (24), 
and remembering that the integration over all space 
has to be performed for the inner space and outer space 
separately, we obtain 

2 Fo 
Asaf (£(u)2u2du 
ky VAn, iua B pas 
nahn ; 2n-+1)); 
+ EA, G jurdu= PPN, (29) 
(27), (28), and (29) now determine completely A, ; and 
a,j; after reduction, and making use of well-known 
integral properties of the Bessel functions, we get 


a os (n+) 
Qa aL 2)}2 Cos T oo) log toe) 
1d 
ie ee ery I) logtvbawt | 
. . « (30a) 
oe (2n-+1) 


ark athl ies es + log ‘eis 


— (gat gag) ecto | 


. . . (308) 
where z= À, jka and y=, kot, 
so that after substitution in (25) we obtain 


spats AZ, | 
Moa 2 ED An, Eb) LD (Ay, kar) Pn (cos $) | 


ky a-o ?O— =) 
(r>a), | 

$ dmi” Te Bn, a 
MaE ae Eo) Pn BIOP, (008 3) ( 
(r<a) 5) 


(31) 
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(31) gives the solution of our problem as a double series 
of Characteristic Functions, the summations to be ex- 
tended over the integer values of n and over all the j values 
every time belonging to a given n. It will be clear that 
(31) must be equal to (1). Further, the orthogonality 
of the functions P, (cos $) requires that for a fixed n 
the summation over the corresponding 7 values must 
every time lead to the coefficients of P, (cos 3) in the 
harmonic series (1), This idea will further be developed 
for the special case |k | =œ in paragraph 7 c. 


Ta. The Harmonic Series for a Perfectly 
Reflecting Sphere. (|ky| =% ). 

In this paragraph we will investigate the limiting 
case for a perfectly reflecting sphere and therefore | ka |=% , 
while basing ourselves on the harmonic series. In the 
next §76 we elaborate the residue series and obtain 
numerical results, showing also the influence of the finite 
height b—a =h of the transmitter above the surface of the 
sphere. In §7c we compare the residue series with the 
development in Characteristic Functions. 

The harmonic series (5) for a=b=r and |+| =% 
becomes 


LiF a i 8). 
Tse 2 ee (2n+1) Nw (cosS)-p,0) + (32) 
where 
1 _ f z aX) (33) 
22Na) | dp aay noS 
ag en (z)} pe 


and we will now examine the numerical structure of the 
coefficient 1/zN(z) as a function of the positive integer 
n for |zi>>1. To this end we can at once make 
use of (22), where, as iz| >> 1, the first term may be 
neglected. We therefore have 


n (2 32 
2N,{z)~—i > G —1 


n ne a) J ® z? 3 
ri gaml 


petir? VED , (34) 
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where 


T z? n 

z < arg V ae 1< 3° 
The square root remains perfectly well defined even for 
z?n? real, when this variable is considered as having an 
infinitesimal positive imaginary part. The approximate 
expression (34) can, when z is not too near n, be further 
considerably simplified, for then the argument of the 
Hankel functions remains so great that the first term of 
their asymptotic development suffices. Under these 
circumstances (34) degenerates to the simple form 


has hl 
NA) Je 
vV 


with the square root in the first quadrant, and this form 
coincides with the first W.K.B. approximation of the 
proper solution of (19) Although (34) was developed for 
n >>I, (35) still gives the right value for n=0, viz. 

-l —i 

2No(z) T : 
Further (35) degenerates for n >>œ>>z to 


(35) 


l œ? 
2N,(z)~ n 
Moreover, for the transition case (n~z) the second 
term of (34) dominates by far, and when n=z exactly ; 
in the latter case this second term becomes indefinite, 
but, while taking the first term of the power series of the 
Hankel functions, it is easy to show that 


{ SA l - } =. 2 IK 1/3) ge~ Zing 3 
zNa z) n=z gue IT(2/3) 


=1:0888... eo 2 ingi, 


The above approximations enable us to determine 
numerically 1/2N,(z) as a function of n in the complex 
plane. We have chosen z=k,a—5°72.10° corresponding 
to waves of 7 meter in the radio case (27a=4.10° em.). 
In fig. 3 we give the result numerically correct. The 
numbers in this figure refer to the values of n—z. It 
appears from this figure that (a) the modulus of 1/zN,,(z) 
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has a maximum near n=z; (b) the density along the 
curve is rather small where the modulus of 1/z N,{z) is 
considerable, whereas nearly all points are lying on the 
two extreme parts of the curve near the two axes; (c) the 
curve has a horizontal tangent at n=z. 

Again, figs. 4 and 5 give the modulus and the argument 
of LzN,,(z) as a function of n—z ; a.o. this clearly demon- 
strates (b) above. 


Fig. 3. 
oso om grien) 


if tz Nal zt 


Coefficient of the harmonic series for |k] =% 
and z=k,a=5-72.108. 


A further numerical consideration shows that for all 
values of n where 1/zN,,(z) has an appreciable value, 7. e., 
where n~z, a good approximation to this function is 
given by the second term of (34) only ; moreover, in this 
range 2?/n2—1 may be replaced by 2(z/n—1), so that the 
final approximation for (34) gives 1/zN,(z) as a function 
of the only variable o=(n—z)/z1”: 


HRZ? 
TEA O C9 


2. 


2N,(2) 


1/35 —2d¢- 
matte tn 
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This shows the universal character of the figures 3, 4, 
and 5, which retain the same form also for other values 
of the wave-length than A=7 meter. 

Returning to the harmonic series (32), the coefficients 
of which have now been completely investigated, it follows 
that, when $ is not near zero, the oscillatory character 
of P, (cos 3) for fixed 4 as a function of n causes the series 


Fig. 4. 
25 
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Modulus of the coefficient of fig. 3 as a function of rs. 


to oscillate very rapidly, so that successive terms nearly 
cancel each other, which makes the summation of (32) 
a very intricate process. We will now turn our attention 
to the much more rapidly convergent residue-series. 


76. The Series of Residues for |k, |=% . 
Consider (7). For |&,| =% the expression N,_1,2 is, 
according to (33), even in n, so that the second term of 
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(7) is exactly zero. We can therefore further use the 
series of residues as given by (15) and for numerical 
applications a knowledge of the zeros of N,_1)/2, as well 


Name at these zeros is necessary. A 


as the values of 
further investigation shows that for |z|>>1, as is the 
case in the radio application, the zeros n, of N,_1,2 as a 
function of n for constant z have a modulus at least of 


Fig. 5. 
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Argument of the coefficient of fig. 3 as a function of n. 


the order |z|. In order therefore to find the proper 
approximation for N,,_1/2(z) we differentiate (16 a) logarith- 
mically, remembering that both |z| >>1 and |r| >>1. 
We thus find 

2 


n z n ain sy du 
naxi ZL {Ten (veri), 
(87) 
Phil. Mag. S. 7. Vol. 24. No. 159. July 1937. M 
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so that the zeros n, of Np- e(z) are given by the trans- 
cendental equations 


itm +4) (°" vie” E E e E acd: 
a1 
(38) 
For real z=k,a (meaning a non-absorbing atmosphere) 


TTT iG 
CCPC RE 
a eneny a 


LTA | tT | 
Yt | 


Poles and characteristic frequencies in the n/z plane for | k, | = 


it follows from (38) that with the approximations used 
all the zeros n, are situated on the curves 


Red 2[" vizet | =o <- - (39) 


in the njz plane, which are shown dotted in fig. 6 as P, 
and P,. The poles n, nearest the real axis and on the 
branch P, with the exception of the first few, are easily 
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found from (38), remembering that |z/n—1| is a small 
quantity. They are given by 
Me=ztre,. ea. > (40) 

where 

Ti=H3r(s +4) Pein, 2 2. (40a) 
In order to calculate the first few roots, which le in the 
immediate neighbourhood of n=z, other authors, e. g., 
Nicholson “2), Laporte), have used approximations with 
more terms. The difficulty is that especially the first 
root, which is of predominant influence on the numerical 
results, lies in the boundary region of the two Debye 
approximations for the Hankel functions; see, e.g., 
Jahnke-Eimde, ‘ Funktionentafeln’ (1933), pp. 204, 207. 
In this way they found, in our notation, 


Nicholson .. 7 )=0°804e'"*, 
Laporte....  7)=0°81 e"%, 
However, it is of interest to note that the expression (36), 


which is valid near n=z, yields as equation for the first 
few zero 


HR (—20)°7} =H { yna) e= bo, 
Remembering further that 


á : e—5/sir 
Hgj(ae-*") = RR a3 {Jarx(a) —J _2)3(«)}, . (41 
and calling «, the zeros of (41), which are all positive, 
and identifying 
22 
32 
we again find for the zeros n, the relation (40), but, 
instead of (40 a), 


(n, —2) 2e Vn — 2T in, 


=}(Bai)t3eim, |... (408) 
which furnishes a better approximation than (40 a). 
Summarizing, the equations (40 a) and (40 6) thus yield 
for the first three roots 


from (40 a): from (40b): 
To= 885e T, T == 0°808 T, 
7y= 2589 eft, TSH DSTT eth, 
T= 5'830 eiT’, E eft, 


M2 
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so that in a relatively simple way we find from (408) a value 
for +, which is very near the accurate values already 
given by Nicholson and Laporte. 
Having thus found the position of the first few poles ns, 
the factor of the corresponding residues, 
oN iz Sie 
(ia), = ee 
may be obtained ee (37) and (40). Substitution of 


the results in (15) yields the final and rapidly converging 
expression 


if Te 143 
II,=4 Be aft ies. y atallar b) ee ene (42) 
[kel = 0 iz 2 s iU D aar, (Ri) Ts 


Taking the special case of the transmitter being placed 
on the surface of the sphere (6=a), (4%) simplifies to 


TT, aap Z Fa), me Ee Bay (48) 
the] = 
where 
3/72 _  pil2r Laren 
F(y)= ws aE : ty (43 a) 


=4/ s 
I= a -a 


D being the distance between transmitter and receiver 
measured along the surface of the sphere. 
The structure of (43), which can be written 


ek D 
Mat Fn), -ooo (AM) 


[esl = 

clearly shows the correction factor F(y) (being a function 
of one variable only and which is due to the curvature), 
which must be added to the solution 2e°/%,D, which 
itself belongs to the plane problem ; the factor 2 repre- 
sents the total reflexion by the infinitely conducting 
second medium; (44) confirms the result first obtained 
by Watson (J. ¢.). 

It is of interest to note that the variable » introduced 
above is closely related with the phase difference ¢ of 
two beams, one going straight to the receiver through the 
sphere, the other going along the surface. In fact for 
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small 5 this phase difference is given by ¢= 7 , and thus 


the solution (44) in this respect shows a great analogy 
with other diffraction problems where also the field near 
the boundary of the optical shadow is expressible as a 
function of the corresponding phase difference. 

From the expressions (43a) and (44) for the vector I the 
normal component of the electric field can simply be 
calculated with the relations given in § 2. 

The result (43 a) is plotted in fig. 7, where the curves 


[E 


Fig. 7. 
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Modulus of the sum of the first few terms of F {y}. 


marked 1,2,3... 6represent approximations of the modulus 
of F(y) as a function of 7 calculated with 1 2, 3... 6 
terms respectively, and where due consideration has been 
given to the phases of the different terms of the series. 

Moreover, referring to (42), where the transmitter was 
not yet lowered to the surface of the earth, this formula 
has been used to calculate the curve labelled h=100 m., 
o=0,A=0-7 metre in fig. 8. As indicated in this figure 
it was constructed for a wave-length of 70 cm. and for 
a height of the transmitter above the earth of h=b—a= 
100 metre. As a comparison we give in the same figure 
the curve marked h—0 for the same wave-length and which 
was calculated from (44). 
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In calculating these curves, as for all following ones, 
the moment of the transmitter has been assumed to be 
such that, when placed on an infinitely conducting plane 
earth, would radiate 1 kw. 


7c. The Characteristic Functions and Values for 
[ka| =o. 
The results of § 6 become much simpler for the special 
case |k,|=0. The |coefficients B, j and therefore 


Fig. 8. 
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Influence of the height of the transmitter field on the 
for 1 kw. “ radiated.” 


also the Characteristic Functions ¢, j then vanish all 
for r<a. Moreover (28), defining the Characteristic 
Values Àn, ;k,, then simplifies to 
d ova TES 
Eewo] a o (0) 
whose roots, according to the definition of {{(z), corre- 
spond to the zeros of the polynomials 


EOE 
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hence, remembering that z{{P(z) satisfies (19), the ex- 
pression (30 a) can be reduced to 


AP (20+ 1) 
mi drka e, Paa 1) — zy 
and the series (31) for b=a to 


tot = =; — = (2n+ 1)P,,(cos i py 


a era er 


vad 
= . (47 
bana j kiltas (47) 


This series corresponds to a point source cos (&,R)/(ik,R), 
representing stationary waves, while the series (32) 
corresponds to a point source e”:"/(ik,R), representing 
diverging waves only. In order to transform (32) to a 
double series analogous to (47) we first write it as 


‘= (2n-+1)P, (cos $). La E (2) | 
a E(k} 


a 
Ma Fa, 
g=kya 


(48) 


Now it may be remarked that the definition of (2), 
as given in § 2, implies that z*+1!e-#{@(z) is a polynomial 
of degree n, and that therefore the last factor of (48) can 
be written as a quotient of two polynomials of degree 
n and n+1 respectively. As the denominator is of 
higher degree than the numerator, this quotient can be 
split up in partial fractions, which, using the differential 
equation (19), can be given the form contained in the 
following series : 


Thov= - X (2n+1)P, (cos 5) 
Kia n—o 


ie ; 1 
oe 2 iw FI) Zi, 3} (ka—an 5)” 
which consists of twice the first part of (47), corresponding 
to the diverging part of the stationary wave represented 
by this latter equation. 
In order to consider (48) in more detail a discussion 


of the roots z,; is necessary, to which we will now 
proceed. 


(49} 
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Like in §7b we therefore have again to consider the 
roots of (37), but this time considered as a function 
of z, n being a positive integer. Therefore we now 
obtain for the equation determining the roots z,, ; instead 
of (38) 


ead i Sa TET le «mjt, G=0,1,...57), (50) 
1 


where, as indicated, 7 goes from 0 to n, because in this way 
we cover all the n+1 roots of (46). For the same reason 
the summation over j in (47) and (49) extends also from 
0 ton. It will be evident that all the roots of (50) are 
situated on the curve in the complex n/z plane given by 


Re f a a2 y So a e a Ü 


u f 
instead of (39). 

This curve is represented in fig. 6 by the drawn sym- 
metrical curve E, on which therefore all the roots z, ; are 
situated approximately, whereas, as shown in §7 6, the 
dotted curves P, and P, contained all the roots n, of the 
equation (38). It is of considerable interest to observe 
that the curves E and P, have in n/z=1 a common point 
and a common tangent making an angle of 7/3 with the 
real axis, as can be easily shown by comparing (39) 
and (51). 

Further away from the point n/z=1 the curve E 
remains finite, whereas P, goes to infinity. In fact, 
by developing (39) for |x/z|>> 1 it can easily be shown 
that the curve P, is asymptotically given by 

p=feltti, 


where p and ¢ are polar coordinates in the n/z plane. 

A further consideration of (51) shows that the curve E 
cuts the imaginary axis in a point A (fig. 6) situated at 
1:509.. =a, say, where « is the real root of 


cotgh (4/143) =/ i+ s 
(see, e. g., Jahnke-Emde, ‘ Funktionentafeln,’ p. 31). It 
can further be shown that for n even there exist 
roots Zn, j=2n, ni2 lying exactly on the imaginary axis 
corresponding to Characteristic Functions containing a 
purely aperiodic time function. The n/z values of these 
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“aperiodic ?” roots converge rapidly towards the value 
given above, as could be verified numerically by deter- 
mining these roots exactly from (46) for small values of 7 ; 
for these small values, however, it is better to consider 
(n-+4)/z instead of n/z, according to (16a). In fact, we 
obtained the following values: 


21/za = 1-5664, 

4h jz, e= 1'526, 

64/zq, 3=1'51Ti, 

8h/zq, a= 1514i, 
Lim n/2,, ye=1-509%. 
i> 


We also calculated the other rootsfor n=1, 2, 3, and 4“), 
and found 


%,0= 0°866—0-5002, Z» o= 1'808—0-702i, 
zı, 1= — 0'866 —0: 500i. 22, 1= 0—1-596i, 
za e= — 1-808 —0°7023. 

2,9 2758—0°843i, |z o= 3°715—0-954i, 
Z= O871—2157i, = | 244 1-752—2-571i, 
z, 2=—0:871—2157 i, 2, 2= 0—2949i, 
Zs, a= —2 758 — 0'843 2. Z, 3=— 1752—2571 i, 
\ (z, 1—3 T15—0 9544. 


The (n+ })/z values of all these roots have also been marked 
with their n value in fig. 6. It is seen that even for these 
lowest values of n» the roots lie already very near the 
universal E, which was derived for n>>1. 

Returning to (50), using a similar procedure as in §7 b, 
we obtain as a general expression for the roots near the 
real axis instead of (40) 


2y,j=n—7m, 2 2... (52) 
where, as in (40 a), 
T= HOm(G+He%e"™, |. (82a) 


Having thus obtained the values of the most important 
roots Z„ j we are now in a position to apply this result 
to transform (49). We therefore invert in (49) the order 
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of summation of n and j and substitute the values of the 
roots 2,,; as given by (52). Retaining the first order 
terms in n, (49) thus becomes 
1°37 =(2n-+1)P,, (cos 3) 
Tig PH z : 
a A nin —ka— (ka) r} 


(53) 


As before, we can transform (53) to an integral over 
n with a path of integration L as in fig. 2c, where now 
the only enclosed pole is 


No=kya + (ka)! Tis 
so that, taking the residue of this single pole, (53) becomes 


oT pno Pui {cos (7 mle 
Tet kya r 7; COS (ngr) 


and, with the approximation (12) for P,, we are led again 
to the series (42), but for b=a. 

This result means that the s’*th term of (42) is the sum 
of the contributions of the Characteristic Functions Øp, s 
where for fixed s the summation extends over all values 
of n from 0 to œ —e. g., the first and dominating term 
of (42) (s=0) is obtained by adding all Characteristic 
Functions „o which are just those which for fixed n 
have the smallest damping. 


8. Strongly Absorbing Sphere 


Uma. V keg? — k? bas } i 
In order to obtain, as in § 7b, the series of residues, but 
now for finite k, a knowledge of the most important 
poles, ż. e., those in the n/z plane nearest the real axis, 
is again necessary. In these poles again N,_12.=0, and 
according to the definition (6) and the approximation 
(37) we have for |n/z| at least of the order unity 


Na-1(2, 2’) 
~— aJ 1. 04% Ny Vel 1%} 


RS ae 8 hog {2'dn(e! ye. (54) 


where 
z=k,a and 2’=k,a. 
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From (54) one can at once deduce the relation which 
must be satisfied by the zeros n : 


imp (™ Vina aa(e+8) 
/,\_% 
V l= 
Z log ey @)} 
dz’ ” 
(s=0, 1, 2, ..-) 


» (55) 


zZ 
z 


where 
0<Re{arcig| J <r. 


Now an approximation, derived in the same way as (35), 
and especially suited when z’—k,a has a great imaginary 
part, becomes 


d r r d + 17H" 
ga 08 bale ~ g l8 2 ae) 


zi ioge meg —i4/1- EE, (66) 
~ da! eae Keg? 


as the zeros n, of N,z_19 appear to lie near z=n. For 
strong absorption [El and |k] are of the same 
order of magnitude, and for the zeros near n/z=1 the 
argument of arc fg in (55) is very small, and hence the 
function arc tg itself, so that as a first approximation 
this term may be neglected in (55). Thus (55) becomes 
nearly equal to (38), 7. e., the zeros near n/z=1 are again 
on the curve E of fig. 6; but a comparison of (55) and 
(38) shows that the zeros are now interspaced between 
the zeros belonging to the case |[k,| =o. This fact 
is especially clear when one considers the zeros in the 
immediate neighbourhood of n/z=1, and which are of 
primary numerical importance. Similarly to (40) and 
(40 a) these zeros are given by 


m=z rt, 2. 2... (87) 

with 
7,=T{3n(s+ Re". . . (57a) 
As in §7b, we again retain in (15) the residues at 


these poles. In order further to find aa for n=n, 
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we make use of (54), and derive from it, Dra (57): 


vV E (—2r AAE a d ogiz p(z") 


oga KA 4t 
: Í ; oN : 
From this expression we derive On for n=n, making 


use of the fact that N,,-1,2=0; thus one obtains three 
terms, of which one dominates for |z|>>1, so that we 
are left with 

ON zfd area 

an ye E log {2'p.(z + | 
Substitution in (15) gives, after reduction, whereby (57) 
is used, ior our final expression : 


ikat 
Taws~ ae kpk AD 


Nx 


Dre iT" 1p? za É; EL (kb) EELT (58) 
ri; Sanh (kya) 
(2=k,a) 
This series shows a great analogy with (42), asy = D/(a?A1*) 
has the same meaning as before. 
If we first limit ourselves to a consideration of the case 
ba, we observe that (58) consists of three factors : 
(1) the first factor 
es jet 
iz tk,D? 
corresponding to the plane problem and for |k,| =% ; 
(2) the second factor 
ikat 
kP (k— k’) D 
representing the effect of the absorption already found 
before for the plane problem; for, with the primary 
field as defined in this paper, the plane problem was in a 
former paper ‘*) already reduced to the exact result 


2k 35 Fah git'Du 1 
Ta tty (oS a aes). 


Rh’ 
at 1 1 
Gam Eat Ey 


and a first approximation for great k,D is obtained from 
this expression by integrating it once by parts and 
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neglecting the contribution due to the upper limit of inte- 
gration, We thus obtain exactly the product of the first 
two factors of (58). 

(3) The third factor of (58) we call G(), and, according 
to the preceding, this factor represents the effect of the 
curvature of the earth. 

Hence, when the transmitter is placed on the surface 
of the sphere, i. e., for b=a, we can write (58) in the elegant 
form 


Tams aes Gq) —=Bne hy Zeta, (59) 
Tavs, a= æ 
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Modulus of the sum of the first few terms of G (n). 


Thus the functions F(q) (of §76) and G(m) (as given by 
(59)) give the correction factors due to the curvature of 
the earth, F(y) being valid for o,=0 and G(y) for 
strong absorption. 

In fig. 9 we have represented by the curves 1, 2,3... 6 
the approximations one finds for |G(y)| if 1, 2...6 
terms in (59) are retained; and thus fig. 9 for Gn) is 
analogous to fig. 7 for F(). As our final result we give 
in fig. 10 the values of the two correction factors |F(y)| 
and |G(m)| applicable to the case o=% and great 
absorption respectively. 

In general, and especially for short waves, not too 
great values of 7 are of primary importance: e.g., for 
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A=7 m., D=50 km., we have y=0-76. Hence practical 
values of 7 are of the order unity, and in this range of 
values of » the difference of |F(y)| and |G()|_ is 
not very great. In practical cases the actual correction 
factor will be found between these two limiting curves ; 
but, e. g., for A=7 m. and sea-water the actual values are 
already very near the curve |G(y)|. This results 
from the fact that then the argument of the are tg function 
(55) is already very small, as the factor 


ka 
(kya PEW ka? —k? 


Fig. 10. 


Correction factors for the curvature of the earth. 


has a small modulus, and the first term of the are tg 
development in (55) suffices for the approximation (59). 
In this way the curves of fig. 11 (for A=7 m.) were 
calculated. In this figure the drawn curves represent 
the electric field of an emitter radiating 1 kw. (for defini- 
tion see §76) for o,=0, for og=10-" and «,=80 (sea- 
water), and for o =10-®, «,=4 (dry soil) respectively. 
The dotted lines represent the corresponding fields one 
would obtain for the same electrical constants of the 
earth, but for a plane surface. This figure shows very 
clearly that for A=7 m. the influence of the finite con- 
ductivity of the soil dominates by far the influence of 
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the curvature of the earth for the distances considered. 
This is a general result for short waves. 

Finally, it was possible with the aid of (58) to caleulate 
the field of an emitter placed at a height h=b—a above 
the ground. This was done for A=7 m., A=100 m., and 
o,=10-15, e,=4. The result is given by the two lower 
curves of fig. 8, where the upper curves gave the result 
for A=0-7m. and o,=«. In both cases no apparent 
effect is seen to occur at the optical horizon of the trans- 
mitter, this being situated at D=35-7 km. A further 
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Influence of the conductivity of the soil on the 
field for 1 kw. “ radiated.” 


investigation showed, however, that for much shorter 
waves (of the order of 1 mm.) a distinct shadow effect 
makes its appearance (cf. Part IL). 
Within a short time we hope to publish Part II. of this 
paper. 
Summary. 


The problem of diffraction of electromagnetic waves 
radiating from a point source outside an absorbing sphere 
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is treated by several methods. In Part I. an exposition 
of the method of Watson is given. This method consists 
in a transformation of the infinite series of spherical 
harmonics for the Hertzian vector into another series, 
which is much more rapidly converging for wave-lengths 
small compared with the dimensions of the sphere. 
Extensive numerical results are given for the case that the 
point source is a transmitter of radio waves, which may 
be situated on the surface of the in general finitely 
conducting sphere as well as at a small height above it. 
The former method is compared with that of development 
in characteristic functions, which is analogous to methods 
used in string problems and wave mechanics. 

In Part II. the original harmonic series for the Hertzian 
vector is split up into an infinite sum of other harmonic 
series. Each of them may be considered as the effect of 
waves (degenerating into rays for extremely high fre- 
quencies) which are once or more times reflected at the 
inner surface of the sphere. These series may be approxi- 
mated by Kelvin’s principle of stationary phase as ex- 
tended by Debye and Brillouin; this method gives appro- 
priate approximations for the field of a radio transmitter 
in the region above its optical horizon, whereas, when the 
point source is at infinity it yields and extends Airy’s well- 
known results for the intensity of the light in a rainbow. 


References. 


(1) Lord Rayleigh, Phil. Mag. (4) xli. pp. 107, 274, 447 (1871); (5) 
xii. p. 81 (1881); (5) xlvii. p. 375 (1899). 
(2) P. Debye, Ann. Physik, xxx. p. 57 (1909). 
(3) Lord Rayleigh, Phil. Trans. Roy. Soc. ceii. A, p. 87 (1904); Sei. 
Pap. F., no. 292. 
(4) A. E. H. Love, Phil. Trans. Roy. Soc. cexv. A, pp. 105-131 (1915). 
(5) G. N. Watson, Proc. Roy. Soc. London, A, xev. p. 83 (1918), and 
A, xev. p. 546 (1919). 
(6) Balth. van der Pol, Phil. Mag. (6), xxxviii. p. 365 (1919). 
(7) O. Laporte, Ann. Physik, lxx. p. 595 (1923). 
(8) T. L. Eckersley, Proc. Roy. Soc. London, A, exxxvi. p. 499 (1932). 
(9) P. S. Epstein, Proc. Nat. Acad. U.S.A. xxi. p. 62 (1935), nr. 1. 
(10) Balth. van der Pol, Physica, ii. p. 843 (1935). 
(11) A. Sommerfeld, Jahresbericht d. deutschen Mathematiker-Vereini- 
gung, xxi. p. 309 (1912). 
(12) J. W. Nicholson, Phil. Mag. xxi. p. 64 (1911). 
(13) Compare: J. J. Thomson, * Recent Researches in Electricity and 
Magnetism, p. 368 etc. 
(14) Balth. van der Pol, Z. Hochfrequenztechn. u. Elekiroakust. xxxvii. 
p. 152 (1931). 


Eindhoven. 
26th October, 1936. 


